Unveiling optical hidden correlations: the illusionist game 
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We suggest and demonstrate an experimental scheme to illustrate and exploit Gaussian correlations 
of quantum states of light. In our setup, termed the optical illusionist game, an illusionist at first 
shows that if two uncorrelated light beams excited in the same Gaussian state are mixed in a beam 
splitter, then no correlations arise between them, as it was not present. On the other hand, when 
correlations with an ancillary state are taken into account, the presence of the beam splitter can be 
unveiled. 



PACS numbers: 03.65.W, 42.50.Dv, 42.50.Xa 

One of the major goals pursued by quantum physi- 
cists is the generation of quantum correlated states for 
realizing protocols with performances beyond any clas- 
sical benchmark PUTS]. These results, besides paving 
the way to new technologies, also provided deeper fun- 
damental insights into the behaviour of quantum sys- 
tems, and the nature of their interactions. In this frame- 
work, Gaussian states of light (e.g., coherent, thermal 
and squeezed light), play a leading role, since they may 
reliably generated, transferred and detected with cur- 
rent quantum optical technology. In fact, using Gaussian 
states, major quantum protocols have been realized, in- 
cluding, e.g., continuous variable quantum teleportation 
[16] and cloning [17]. Linear optical elements may be in- 
deed used to create and manipulate correlations of Gaus- 
sian states [TB1 - E3"] . and their use has been extensively 
analyzed in the recent years piH?f] , in order to optimize 
the generation of correlations, and for their preservation 
during the propagation through noisy channels. Notice 
that even if Gaussian entanglement is the most impor- 
tant resource for continuous variable quantum technolo- 
gies, more general quantum correlations, as the Gaussian 
discord [351 [55] , may represent a tool for enhanced pro- 
tocols as well [30l [31] . 

In this Letter, we address and exploit the invariance 
of identical (uncorrelated) Gaussian states mixed at a 
beam splitter (BS). In particular, we demonstrate that 
they are left locally unchanged, as it happens when no 
BS is present, but nevertheless one can discriminate be- 
tween the presence or the absence of a BS, uppon ex- 
ploiting the correlations with a third ancillary mode. 
In order to present the phenomenon to a general au- 
dience, we illustrate the protocol in terms of a game, 
which we term optical illusionist game, and which pro- 
ceed as follows: An illusionist displays on the stage an 
amazing game involving two light beams, a beam split- 
ter (BSi) and two photodetectors arranged to measure 
the correlations between the two beams. He, or she, 
specifies that only Gaussian states of light will be used, 
with the additional requirement of having zero (coher- 




FIG. 1: The optical illusionist game, (a) Two uncorrelated 
light beams excited in the same Gaussian state are crossed: no 
correlations arise, (b) A beam splitter (BSi) let the two beams 
interfere: still no correlations are detected, (c) Nevertheless, 
the illusionist can always say whether the BSi is hidden under 
the hat or not: the trick is behind the curtain, where the 
illusionist has access also to beam 3, which is generated in a 
correlated state with beam 2. Dk, k = 1,2,3, are detectors 
to measure the correlations. 



ent) field amplitude. The illusionist starts the perfor- 
mance by crossing two uncorrelated beams of light ex- 
cited in the same Gaussian state and showing that no 
correlations arise between them (Fig. [I^,). Then, he in- 
serts BSi and shows that the detectors continue to re- 
veal zero correlations (Fig. [T]d). This effect is a conse- 
quence of the linear interaction of modes performed by 
the BSi and the Gaussian nature of the states [211 [25] . 
In fact, the most general single- mode Gaussian state can 
be written as g = D(a)S(£)v th (N)S*{£)Di(a), where 
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vth(N) = E"=o( iV ) I V(l + N) n+1 \n)(n\ is a thermal 
equilibrium state with N average number of quanta, 
S(r) =exp{![£(at) 2 -£*a 2 ]} andD(a) = exp[atf -a*a] 
are the squeezing and the displacement operators, respec- 
tively, a being the annihilation operator. The state g is 
fully characterized by its 2 x 2 covariance matrix (CM) 
<x, [cr] hk = \{R h R k + R k R h ) - (R h ){R k ), k = 1,2, and 

first moment vector X = (.R T ) = v / 2(3fJe[a], 3m[a]), 
with (A) = Tr[Ag], where R T = [R 1 ,R 2 ) = (q,p) 
and q = (a + a))/\/2 and p = (d — a 1 * ) / (iy/2) are the 
position- and momentum-like operators. In the follow- 
ing, according to the requirements of the performer, we 
— T 

set a = 0, i.e., X — 0. Now, when two uncorrelated, 
single-mode Gaussian states g k with CMs <r fc , k = 1,2, 
interact through a BSi with transmissivity t, the initial 
4x4 CM S = <ti © cr 2 of the two- mode state pi ® p 2 
transforms as I32TI341: 



correlations between modes 2 and 3. Since the BSi acts 
only on modes 1 and 2, we have: 
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where Si = tct\ + (1 — r)er 2 , S 2 = t<t 2 + (1 — t)ct\ and 
S12 = t(1 — t)(<t 2 — <7i). Note that Si 2 7^ denotes the 
presence of correlation between the outgoing modes. It is 
now clear that if the input states are in the same initial 
state, namely cr k = er, then S = S*- "*- 1 , i.e. the output 
beams are left uncorrelated (Si 2 = 0), independently on 
the trasmissivity of the beam splitter. 

In the second part of the game, the public is requested 
to insert at will the BSi hiding it under a hat (Fig. [TJ;). 
The illusionist argues that he is able to say whether it 
is present or not and he succeeds, despite the fact that 
the interaction through the BSi appears to not affect the 
output states. 

How is this possible? The trick stays in the use of a 
third ancillary mode, which is correlated with one of the 
two input modes by the initial preparation, and allows 
the illustionist to discriminate between the two cases. In 
fact, although the public sees only modes 1 and 2 on the 
stage excited in the state gi = g 2 = g, mode 2 is actually 
correlated with mode 3 (kept by the illusionist behind the 
curtain, Fig. [IJi), i.e., g 2 = Tr[p 23 ] with g 23 ^ g 2 ® Q3- 
On the one hand, modes 1 and 2 are left unchanged and 
uncorrelated after the interaction at the BS. On the other 
hand, because of the interaction, part of the correlations 
shared between modes 2 and 3 are now shared between 
modes 1 and 3. This can be seen by looking at the evolved 
CM of the whole state of the three modes in the presence 
of the BSi. The 6x6 CM of the initial state gi 23 — 
Qi ® Q23 reads: 
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where cr k is the 2x2 single-mode CM of mode k = 1, 2, 3, 
a 1 = a 2 = cr and the matrix <5 2 3 7^ contains the 
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The comparison between Eq. ^ and Eq. ^ shows that 
while the states of modes 1 and 2 are (locally) left un- 
changed and, in turn, uncorrelated, both modes 1 and 2 
are now correlated with mode 3 (again, the presence of 
non-zero off-diagonal blocks denotes the presence of cor- 
relations between the corresponding modes). Further- 
more, the degree of correlations between the modes 2 
and 3 is decreased (8 23 — > %Jt 6 23 ) for the benefit of the 
birth of correlations between the previously uncorrelated 
modes 1 and 3 (0 — » y/i — r <5 23 ). We conclude that BSi 
is present, since mode 1 becomes correlated with mode 3 
at the expenses of the initial correlations between mode 1 
and 2. Nevertheless, this happens in such a way that no 
overall correlations arise between the interacting modes, 
but the illusionist has also access to the third mode and, 
thus, can discover the presence or not of the BSi: he 
exploits "hidden correlations" . Remarkably, the correla- 
tions are not only hidden "behind the curtain" , to use the 
terminology of our fancy game: at the BSi the two beams 
are actually mixed (part of the photons from each beam 
is reflected and transmitted according to the interference 
rules of quantum physics) and, shot by shot, correlation 
may arise due to the quantized nature of the light, but, 
overall, this effect is "washed out" by the very statistical 
nature of the Gaussian states and the linear interaction 
itself. Thus, the ancillary mode allows retrieving this 
particular information. 

For the sake of simplicity, we have described the illu- 
sionist game by addressing the use of single-mode beams. 
Notice, however, that the same results hold also in the 
presence of multimode Gaussian states. This is essen- 
tially due to the tensor product nature of the multimode 
state and the pairwise interaction, which ensures that 
each mode interferes with the corresponding one. 

We have experimentally realized the protocol by ad- 
dressing the most simple Gaussian states, i.e., thermal 
beams (Fig. [2]). Two speckled, i.e., spatial multimode 
and single temporal mode pseudo-thermal beams, a and 
b, were generated by scattering two 1 ns laser pulses 
(A = 532 nm and 12.4 Hz repetition rate) on two inde- 
pendent rotating ground glasses (Rl and R2) [36]. While 
beam b is directly sent to BSi as beam 1 (Fig. [2^,), beam a 
is divided into the two correlated beams 2 and 3, needed 
for the trick, by means of the balanced beam splitter 
BS 2 3. A system of half waveplates (A/2) and polarizers 
(P) compensated unbalancing in intensity of the modes 
and selected only modes with the same polarization, in 
order to realize the theoretical condition g\ = g 2 . The 
lens LI (focal length 20 cm) imaged the beams on a EM- 
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FIG. 2: Sketch of the experimental set-up. Implementation of 
the illusionist game with the BSi (a) and without it (b). In the 
right panels of the figures (a) and (b) the thermal beams and 
Ak are reported: notice that the presence of BSi reduces the 
correlations between areas A2 and A3 (a), which are apparent 
when it is not present (b). 



CCD array (Andor Luca, sensitive area 658 x 498 pixels, 
squared pixel dimension 10 fim) with 0.16 magnification 
of the beams. Each mode of the beams corresponded to a 
point on the EMCCD array. For each beam k — 1, 2, 3 we 
selected an area Ak with 40 x 40 pixels collecting M w 41 
modes. Each mode collected in the area A/, interferes 
with its correlated one in the corresponding other areas. 

The pseudo-thermal beams presented the typical 
speckle structure; we put the BSi in their far field zone in 
order to associate each mode of the multimode thermal 
beam to a point on the BSi, i.e., on the detection plane: 
this setup allowed identifying properly each pair of single 
modes and interacting them at the BSi. 

To evaluate the intensity correlations between the 
modes, which are due to the presence of non-zero off 
diagonal blocks of the CM, each beam k = 1,2,3, was 
registered by an area Ak of the EMCCD array (right 
panels of Fig. [2]), corresponding to the detector Dk of 
Fig. [l] Due to the presence of many spatial modes per 
area Ak, we actually measure the corresponding inten- 

rh) = E 




FIG. 3: Plot of the correlation ci,2 for different experimental 
configurations (dots), together with the confidence intervals 
at 99 % (shaded regions), as function of the number of frames 



used to evaluate it. (black) and cf^ (orange) correspond 
to the correlation of the beams 1 and 2 outgoing BSi when 
only mode a or b, respectively, is switched on. c 1 B 2 S ' ) (blue) 

and Cj 2 (red) refer to the correlation of modes 1 and 2 with 
and without the BSi, respectively, when both a and b are 
switched on. As expected, the width of the confidence inter- 
vals decrease as increasing Nf Tames increases. The insets of the 
figure show the probability density function of the correlation 
coefficients for N{ Tanlc — 10, 30 and 50. 



Ch,k = 1, one has maximum correlation; if Ch,k = 0, the 
beams are uncorrelated. It's remarkable to note that if 
all the spatial modes belonging to a beam have the same 
intensity the correlation coefficient c% k is independent of 
the number of modes M, i.e.., Ch,k is equal to the one 
calculated for the single mode. To ensure proper mode 
matching at BSi and perfect superposition of the corre- 
lated areas A\ and A2, we alternatively stopped beams 
a or 6 and we measured and maximized the correlation 
and c^2i respectively, between the beams 1 and 2 
outgoing the BSi. We then investigated the two config- 
urations proposed by the illusionist: with the BSi and 
without it. Our results fully demonstrate the theoretical 
protocol. Fig. [3] shows the correlation coefficient c 12 for 
different configurations and as a function of the number 
of frames. With the BSi inserted, the correlation falls 
from almost 1 in the presence of only one of the two 



Ch,k 



<Kn,k°™,k) for each frame j taken with beams a and b (c[ a j and c[ > 2 ), to almost zero when both 

beams are present, thus efficiently hiding the correlation 
(q 2 ). A very small residual correlation persists due 
to experimental imperfections, such as the non perfect 
superpositions of the optical modes and the unbalanc- 
ing of the BSi. The question that naturally arises is 
whether this residual correlation is enough to discrimi- 
nate the presence from the absence of the BSi. The an- 
swer is "no" where a relative small number of frames is 
used to estimate c\^- As showed in the insets of Fig. [3j 
the large overlap between the probability density func- 
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the EMCCD, where a m .k is the field operator of the Til- 
th mode impinging on Ak and M is the total number of 
collected modes. The correlation between the beams h 
and k (h,k — 1,2,3) is estimated by using the second 
order correlation coefficient c^kt namely: 



(jfcjft)fr — (Ih)tr(Ik)h 

A lY {I h )A h {h) 
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where (F) b = (N iia , mc ) 1 Z),=T m ° is the average 
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FIG. 4: (a) Plot of ci,s (dots), together with their con- 
fidence intervals at 99% (dashed lines and shaded regions), 
as functions of JVf ramc . c^-p (blue) and c[ n ^ (red) refer to 
correlations between beams 1 and 3 with and without BSi, 
respectively. (b)Plot of C23 (blue) and c^' (red), referred 
to correlations between beams 2 and 3 with and without BSi, 
respectively, (c) Log-linear plot of the error probability Pjrr ^ 
of the discrimination presence/absence of BSi as a function 
of iV fr amc 

obtained from Chk- 



-Nframe = 10,30,50) does not allow to discriminate be- 
tween the two cases with sufficiently small error probabil- 
ity and the "hidden correlations" are needed. In Fig. [4^ 
we report with and without BSi: the correlations 
between beams 1 and 3, which are initially uncorrelated, 
are evident in the presence of BSi. In the absence of BSi, 
beams 1 and 3 are uncorrelated, hence c^ 1 ^ « 0; inserting 



BSi, correlations appear and c\ 
one has c 
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0.5. Similar results 
s 1 without BSi and 



c 23 ' ~ 0-5 w ith BSi. However, to perform the trick, the 
illusionist should be able to discriminate the two cases 
with low error probability and, thus, he should choose 
a suitable Nf lamc . In Fig. [Ib we plotted the probability 

(1 2) 

of error P er r as expected by the public, i.e., addressing 

( 1 3) ( 2 3) 

ci 2 , and that of the illusionist, 

Perr or Pcrr , i.e., ad- 
dressing ci t 3 or 02,3, respectively. Also in the presence 
of experimental imperfection, for A^ rame ss 50, he is in- 
deed able to discriminate between the two configurations, 
since P c ( r V 3) w 1% ( P ( r 2 r 3) « 1%) and P^ 2) » 35%: we 
conclude that the illusionist succeeded in his game. 

In summary, we have suggested and demonstrated a 
novel scheme to illustrate and exploit the counterintuitive 
correlations of quantized light. In our illusionist game 
two identical Gaussian beams of light are addressed to 
the BSi and do not acquire any correlation revealing the 
action of the BS itself. On the other hand, the presence of 
the BS can be identified by using the pre-existing correla- 
tions of one of the beams with a third one. Our protocol 



represents a paradigmatic example of hidden correlations 
of quantized optical beams, and paves the way for future 
applications, e.g., innovative interferometric schemes and 
realization of high averaged photon number entangled 
states [3"fj]. 
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